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Abstract 

I present a theoretical model of a quantum statistical ensemble for which, unlike in conventional 
physics, the total number of particles is extremely small. The thermodynamical quantities are 
calculated by taking a small N by virtue of the orthodicity of canonical ensemble. The finite 
quantum grand partition function of a Fermi-Dirac system is calculated. The model is applied to 
a quantum dot coupled with a small two dimensional electron system. Such system consists of an 
alternatively single and double occupied electron system confined in a quantum dot, which exhanges 
one electron with a small N two dimensional electron reservoir. The analytic determination of the 
temperature of a (1 2) electron system and the role of ergodicity are discussed. The generalized 
temperature expression in the small N regime recovers the usual temperature expression by taking 
the limit of N — > oo of the electron bath. 
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I. INTRODUCTION 



A deep connection exists between the nature of the elementary objects described by quan- 
tum mechanics, and the emerging properties of thermodynamical quantities. Bohr clarified 
that the complementarity principle should apply to energy and temperature measurements. 
Indeed, the determination of the former is incompatible with the knowledge about the sec- 
ond, when they refer to elementary objects lfl. After the creation of solid state quantum 
dots [2], it became possible to explore confined fermionic systems constituted of few elec- 
trons p(, down to a single localized electron Such systems are electrically probed by 
means of accurate charge sensing capable to determine current fluctuations corresponding to 
a variation of charge much below the charge unit js-ll]. The experimental determination of 
the electron temperature in a nanostructure is a difficult task [121 ]. but it becomes a principle 
issue when the system consists of a (1 -H- 2) electron system in particle exchange with a few 
electron bath. Even if the definition of temperature T = ff^j) in terms of the energy 
U and the entropy S holds independently from the size of the ensemble 131 ] . the systems 
here considered are far from the conventional assumptions to derive thermodynamics from 
statistical physics. Apparently, it should not be possible to associate a temperature to such 
small systems. The reduction of the size of the system down to one electron in thermal 
and particle exchange with a finite electron system (small number of electrons N) implies a 
twofold change of perspective to determine the thermodynamical quantities. The first major 
change is the shift from space to time ensembles. Statistics can be recovered on a few particle 
open system only by considering the average of measureable quantities in the time domain. 
The second change consists in the generalization of thermodynamics by removing the limit 
of large N, which implies that the terms at the order 1/N are relevant in the determination 
of the physical quantities. The equivalence under the first change of perpective is granted 
by the ergodicity, while the second change is granted by the orthodicity. In the Section II, 
I define the finite quantum grand partition ensemble in the limit of small N. Next, I deter- 
mine the occupation probability of an electronic system capable to contain either one or two 
electrons. The state equation is derived and it surprisingly depends on the heat capacity per 
area unit ca- In Section III, the determination of the generalized temperature expression 
and its dependence from experimentally measurable quantities is described. Such evaluation 
is done by considering the realistic experimental condition of a single electron quantum dot 
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coupled with a small electron reservoir. The results are summarized and briefly discussed 
in the Conclusions. 



II. FINITE QUANTUM GRAND CANONICAL ENSEMBLE OF A (1 < — ► 2) 
FERMIONIC SYSTEM 

The grand canonical ensemble consists of a large open ensemble made of identical systems, 
which is in thermal equilibrium with a reservoir at a given temperature. The ensemble 
under investigation and the thermal reservoir may exchange energy and particles. In this 
section the statistical physics of a similar system is investigated. The model is derived for a 
generic system of fermions and it applies to electrons in a quantum dot. Contrarily to the 
grand canonical ensemble, here the open system under investigation may contain only 1 or 
2 fermions. The reservoir is a low dimensional system of few fermions. According to the 
physics and technology of semiconductor nanodevices, to which the present study can be 



applied, the reservoir is well approximated with a two dimensional electron system (2D 



DSj, 
141. 



with a negligible extension in the direction perpendicular to the plane of the electrons [?]• 
To grant the uniformity of the notation, the island capable to capture one extra electron 
is exemplified with a trap spatially extended along two dimensions, without any loss of 
generality. Such assumption is realistic, being the confinement of the electron wavef unction 
symmetric in the x — y plane both in defects at the Si / Si02 interface and in a quantum 
dot which is either litographically or split gate defined in a heterostructure 14j. Therefore, 
all the definitions in the following refer to surfaces Sj instead of volumes Vi where i = 1 
indicates the quantum dot and i = 2 a thermal bath made of N 2 electrons. 

A system is orthodic if the averages of the physical quantities on a given distribution 
recover the laws of thermodynamics. The canonical ensemble is natively orthodic so such 
properties is granted for the finite N grand canonical ensemble by construction. 

Said H(N) the Hamiltonian of the system constitued by N electrons distributed between 
the island and the 2DEG, it can be separated by H = Hi(Ni) + H 2 (N 2 ) where Ni = 1, 2. 

Consequently one defines the finite quantum grand partition function for identical parti- 
cles: 
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Z(E,N,T) = ^(S 1 ,7V 1 ,T)Z(E 2 ,7V 2 ,T) = (1) 

JVi = l 

= j^tr (e-^) ir (e^ 2 ) (2) 

AT 1= 1 

where iV = Ni + N 2 and E = S x + E 2 . It is introduced the function 

P(7Vl) "Z(S 1 + S 2) iV 1 + iV 2 ,T) e (3) 
where /3 = (fcT) -1 which by definition satisfies: 

E^(p(iV 1 )) = l (4) 

JVi=l 

Since tr (p) is the canonical partition function times Z(E 2 , N 2 ,T)/Z(E, N,T), in the 
following such ratio is calculated by means of the Helmholtz potential \1/ = —(3~ l logZ, 
which gives: 

Z(T, 2 ,N 2 ,T) _ / 35r( Sl+S2i iv 1 +Ar 2 ,T)- / 9*(E 2 ,Ar 2 ,T) /r\ 

Z(E,N,T) K) 
In order to explicitly evaluate the internal energy U is now calculated, since \1/ = 
U-TS = U + T§; where S = - (ff ). Temperature T is defined as usual from T = (|g )~\ 
For a Fermi 2DES of M electrons in a surface with area A, in the limit of small temperature 
T, it holds: 

U(M, A) = A ■ u = A Eg(E)dE + ^ (k B T) 2 {pig' (ft) + <?(//)] j (6) 

M = A-n = A^£ F g(E)dE + ^ (k B T) 2 [«/(/.)]} (7) 

where the density of states per surface unit is g(E)dE = -^j^dE for a two dimensional 
system with d — 2 and //p is the Fermi energy. Since g'(E) = at d = 2, 

f/(M,A) = |A/i| + ^(A; B T) 2 (8) 

M = gAfip (9) 
The inversion of Equation 9 gives the chemical potential 

MM,i) = ^ (10) 



so the internal energy U is equivalently written 

1 M 2 it 2 

U(M) = -— + -gA(k B T) 2 (11) 

The above relations are valid for M = iV 2 electrons with A = E 2 and for M — N± + N2 
with A = Ei + E 2 , while U(Ni) is treated separately because of the dependence of the 
nature of confinement when Ni = 1,2. The polynomial shape of U(M) = au,A + oaT 2 , 
where o,m,a = an< ^ = ] ^gA{k B ) 2 as a function of the temperature T implies that 
*(M,A) =a M ,A-b A T 2 so 

1 M 2 it 2 

nM,A) = -—--gA(k B T) 2 (12) 
It is useful to calculate the heat capacity per area unit at constant surface 



and the pressure 



(dV(M,A)\ IM 



2 7T 2 



P ( M '^ = -(^^J T =2^ + Y^ 2 < 14 > 

Since the island which traps the electron has a very small spatial extension, it holds the 
approximation E = E2 which simplifies the analysis. It is now possible to evaluate the 
Helmholtz energy variation between JVi + iV 2 and iV 2 electrons: 

A* = + 7V 2 , E x + E 2 ) - *(7V 2 , E 2 ) = (15) 

■^S-a-^^^^*^"- (16) 

7V lM (7V 2 ,E 2 )(l + ^)-l CA EiT (18) 



A\I> goes in the Eq. 5 and gives 

Z{T. 2 ,N-N U T) = M Ar 2 )iv 1 (i + ^)-^ 

Z(E,7V,T) 1 J 

The /mite quantum grand canonical ensemble can be therefore defined by the pair (£,p) 
with £ = T and: 
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where z = e^ 1 ^ 2 **, while T indicates all the states of the system over which the summations 



are performed and it consists of the discrete quantum analogue of the Gibbs r-space. [151] 
The probability of occupation of the subsystem 1 with L fermions is therefore 



trp(L) z W,2N 2)tr ( e -^ l( p- 



z L(l+L/2N 2 ) tr ( e -pH x {L) 



(22) 



Zqg 
where 

Z QG = £ ^(i+m/wifr ^-mW} (23) 

iVi=l 

The present section is concluded by expressing the state equation and the relationship 
between the temperature T and the occupation statistics. Since J2"n 1= i t r (p(-^i)) = 1> h 
holds: 

Z QG -e~^F = l (24) 
or equivalently, if one considers the logarithm of both sides: 

logZ QG = (25) 
The generalized temperature of a (1 •H- 2) system is therefore given by inverting 

p(l) _ /^tr^^W) ^ 

^ " /( 1+ *)tr(e-^( 2 )) 
In the next section such ratio is explicitly evaluated by considering the realistic value of 

a possible solid state quantum device. 



III. SINGLE ELECTRON TEMPERATURE IN A QUANTUM DOT WITH A 
SMALL N ELECTRONS RESERVOIR 

In this section the physical parameters involved in the electron occupation probability of 
a realistic quantum dot and the consequent experimental determination of the generalized 
temperature are discussed. The study of the electron occupation of a quantum dot can be 
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realized by measuring its charge state by the current in a channel electrically coupled to 
the electron charges confined in the dot [si, If]]. In the case of a natural quantum dot like 
a donor or a lattice point defect close to the Si/SiOo interface, the channel is provided 

n n 

by the two dimensional gas formed at the interface by applying a gate voltage |9|, II If . 
In the case of lithographically defined quantum dots, a current flows in the proximity of 
the island which confines the localized charges Let's consider the simplest system as 
possible, like a point defect close to the Si/Si02 interface. Typically the 2DES is confined 
in a (50 — 300) x (50 — 300)nm 2 , while the electron wavefunction is spread along few nm 
in the direction perpendicular to the 2DES [181 ] . The point defect can only accept one 
extra electron. When the defect is paramagnetic (it switches from jVj = 1 to N% = 2 
and viceversa), the first electron fills the ground state of a hydrogen-like shell, at energy 
E(l) = Et- Indeed, the high extraction energy of the unpaired electron makes impossible to 
achieve the ionization of the first localized electron, unless a metal gate is used to manually 
modify the charge state from Nx = 1 to N± — 0, differently from the case here considered. 
Because of the spin degeneracy, a second electron can be captured at the same energy level 
Et to constitute a singlet state ^= (|t) \i) — \l) |T))? with two extra contributions due to the 
Coulomb charging energy AEq and to the lattice relaxation AE^. While the origin of the 
first contribution is straightforward, the second requires a short discussion. The presence of 



the second electron involves indeed a rearragement of the lattic e |17l| . At the low temperature 
here considered the 2DES is weakly coupled with the crystal 19| . Phonons are involved in 
the emission and capture of one electron from the defect and the relaxation of the crystal 
implies a change of energy of Shr^ where Shr is the Huang-Rhys factor and u is the 



average phonon frequency in the configuration coordinate picture [20j]. In the following we 
consider AEl = ShrFioj the energy gain of the lattice when an electron is captured. In a 
natural quantum dots it is not possible to capture a third electron because the energy of the 
system would exceed the conduction band edge energy. 
The total energy of the two electrons is therefore 

E{2) = 2E T + AE C + AE L (27) 

Experimentally, the time resolved trace of the current is analyzed by means of fast dig- 
itizers and they typically appear as a random telegraph signal for both the quantum dots 



161 ] and the point defects 



lOt ] . A typical switching current trace is shown in Figure 1. The 
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capture time r(l) and the emission time r(2) are obtained from the average of thousend of 
switching events and they generally obey to a Poissonian statistics. Deviations from such 
statistics are quantified by a non unit Fano factor and they reflect quantum coherence typ- 



ical of pecularly coupled quantum dots 



211 ] . Their ratio is governed by the ratio between 



the occupation probability between the two occupation probabilities: 

r(l) p(l) 



(28) 



r(2) p(2) 

It is consequently suitable to calculate the relative occupation probability for the two 
states with 1 or 2 electrons, since their ratio can be experimentally accessed. For the 
particular case here discussed, when the subsystem 1 is a paramagnetic point defect 

m = 2e M i+ ^) e ~ PET = (29) 

p (2) e -P2E T +AE c +AE L 



i(e t +ae l +ae c -ij,(n 2 )(i+ 1 



The factor 2 takes into account the spin degeneracy of the ground state. The values of 
p(l) and p(2) are experimentally obtained as the average occupation time in the states 1 
and 2 respectively. Such identification is possible by virtue of the ergodic hypothesis. The 
complete experimental determination of the parameters involved in the Eq. 30 provides the 
generalized temperature T shared by the electrons in the quantum dot and the electrons in 
the 2DES. The generalized temperature of the electron(s) localized in the dots in thermal 
equilibrium with the small N 2 electron bath may consequently be defined: 

T = tfln^Q (Br + AE L + AE C - fi(N 2 ) (l + JL)) (31) 

It is remarkable that such definition of temperature of a time ensemble of few electrons 
concides with the usual one for a space ensable just by taking N 2 — » 00. For this reason it 
can be considered a meaningful extension of the temperature definition for a small quantum 
system of electrons of the kind treated in the present paper. Such result can be easily 
adapted to the case of the confinement induced in a artificial quantum dot, which does not 
require the lattice relaxation term. 
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IV. CONCLUSION 



The finite quantum grand canonical ensemble (£, p) has been defined for a fermionic 
systems constituted by a cell capable to confine either 1 or 2 electrons, and a thermal bath 
of finite and small size made of a two dimensional system of few electrons. The state equation 
is governed by the heat capacity per surface unit ca instead of the pressure. Averaging on 
a time ensemble substitutes the average on space ensembles for the determination of the 
thermodynamical quantities. The ensemble has been applied to a quantum dot constituted 
by a natural point defect at the Si/Si0 2 interface and it holds for a general (1 -H- 2) system. 
The ratio of the characteristic times monitored by the current two-state fluctuations is given 
by the ratio between the occupation probability calculated with the presented approach. 
Therefore, the generalized temperature of such a small time ensemble can be defined and 
extracted from the ratio between the average characteristic times of the two current states. 
Such definition of temperature returns the usual temperature by taking the limit of N — > oo 
of the electron bath. 
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V. FIGURE CAPTIONS 



Fig.l The change of the charge occupation of both an artificial quantum dot and a 
natural point defect is monitored by recording the sudden changes of an ultra 
weak channel current electrically coupled with it. The capture and the emission 
phenomena are generally governed by a Poisson process which determines the 
average occupation time. Ni is the number of electrons in the dot. In this 
example the high current state is associated to the N± — 1 occupation. The 
two states could also be reversed (Ni = 2 would refer to the high current state) 
depending on the microscopic nature of the electrostatic coupling of the island 
with the two dimensional system. 
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